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ABSTRACT 

We study C/O white dwarfs with masses of 1 .0 to 1 .4 Mq accreting solar-composition material 
at very high accretion rates. We address the secular changes in the WDs, and in particular, the 
question whether accretion and the thermonuclear runaways result is net accretion or erosion. 
The present calculation is unique in that it follows a large number of cycles, thus revealing 
the secular evolution of the WD system. 

We find that counter to previous studies, accretion does not give rise to steady state burn- 
ing. Instead, it produces cyclic thermonuclear runaways of two types. During most of the 
evolution, many small cycles of hydrogen ignition and burning build a helium layer over the 
surface of the white dwarf. This He layer gradually thickens and progressively becomes more 
degenerate. Once a sufficient amount of He has accumulated, several very large helium burn- 
ing flashes take place and expel the accreted envelope, leaving no net mass accumulation. 

The results imply that such a system will not undergo an accretion induced collapse, nor 
will it lead to a SN Type la, unless a major new physical process is found. 

Key words: novae, cataclysmic variables — supernovae; general 



. 1 INTRODUCTION 

The prevailing scenarios leading to Type la SNe can be divided 
, into two classes - those of singly (SD) and those of doubly de- 

■ generate (DD) systems. In the first, a degenerate WD accretes 
] from its binary companion, and accumulates sufficient mass to ap- 
. proach the Chandrasekhar limit, where CO detonation or the col- 

■ lapse to a neutron star takes place (e.g., Hillebrandt & Niemeyer 
, 2000; Podsiadlowski et al. 2008). In the second class, la's occur 

■ from the merger of a binary WD system (e.g., Webbink 1984). 
' Since no model within the above classes of progenitors comes with- 
, out significant caveats, there is still no consensus model for Type 
• la SNe. 

, When considering the merger of WDs, one major caveat 
is the merger rate. The theoretical and observed rates of super- 
Chandrasekhar mergers is typically an order of magnitude smaller 
than the observed rate of la's, unless carbon detonation can take 
place in sub-Chandrasekhar systems as well (Van Kerkwijk et al. 
2010; Badenes & Maoz 2012). Another problem is that merger sce- 
narios would tend to produce explosions which are more heteroge- 
neous than observed la's, both in terms of light curves and ejected 
mix of elements. 

Systems where a WD accretes from a non-degenerate com- 
panion, on the other hand, are subject to other theoretical and ob- 
servational limitations. Perhaps the most important one is that of 
the accretion rate. In their review, Hillebrandt & Niemeyer (2000) 
pointed out that although SD is the favoured progenitor model, the 
major problem has always been that nearly all possible accretion 
rates can be ruled out by strong arguments. 



When hydrogen is accreted at sufficiently low rates, it can cool 
and become degenerate. When a sufficient amount of gas accumu- 
lates, the hydrogen ignites to produce an unstable thermal flash 
(Schwarzschild & Harm 1965; Rakavy & Shaviv 1968). Since the 
nuclear luminosity may reach extremely high values, of the order 
of lO'^Lo for a relatively long period of time, the flash can dynam- 
ically eject the accreted mass. In fact, this mechanism naturally ex- 
plains classical novae eruptions (Iben & Tutukov 1984). Because 
lower accretion rates allow more mass to accumulate, lower accre- 
tion rates produce stronger eruptions. Moreover, the observations 
indicate that the ejecta contain as a rule, more heavy elements than 
in the accreted matter. It is obvious that if He burning and beyond 
does not take place, then these heavy element originate form the 
underlying WD, which is eroded in this process. 

Prialnik & Kovetz (1995) and Yaron et al. (2005) found that 
all WD's accreting at rates m ;£ lO^'Mo/yr erode in mass. At some- 
what higher accretion rates, there is net accumulation, but with a 
low efficiency since most of the mass is still ejected. 

One way to avoid the hydrogen flashes associated with the 
low accretion rate, and the consequent mass loss, is to accrete pure 
helium, which may avoid ignition (Iben & Tutukov 1991, 1994). 
However, population syntheses suggest that such progenitor sys- 
tems, which appear as AM CVn stars can contribute at most, about 
10-- of the SN la rate (Solheim & Yungelson 2005). 

At the opposite limit, of very high accretion rates, 
Iben & Tutukov (1984) argued that the WD does not experience 
hydrogen shell flashes, but instead, its envelope expands to a radius 
R « 1050(M/Mo - 0.5f-^^Ro, which for a massive WD is much 
larger than the orbital separation of typical cataclysmic variables. 
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Once a common envelope is formed, heavy mass loss should then 
prevent the mass of the C/O WD from reaching the Chandrasekhar 
limit. 

In between the above ranges, are accretion rates for which 
the released gravitational binding energy is close to the Edding- 
ton limit, around 10"' to IQ^^ M^lyr (depending on the WD mass). 
It was suggested that the high accretion should lead to a quiet hy- 
drostatic burning of H and He (Nomoto 1982). 

Various attempts to calculate the effect of high accretion rate 
were carried out under different approximations (e.g., Nomoto 
1982, using the technique developed by Nomoto & Sugimoto 
1977). In particular, the assumption of steady state ignores the 
time dependent evolution of the accretion system. However, this as- 
sumption simplifies the calculation because time dependent models 
have two major obstacles. First, the latter require very fine mass 
shells, and second, the full evolution requires following a large 
number of small flashes. 

Fujimoto (1982) investigated the thermal properties of hydro- 
gen shell burning on accreting white dwarfs. Assuming hydrostatic 
equilibrium (no dynamic effects were allowed), Fujimoto found 
that the hydrogen bums steadily on the surface of the WD. The 
gravitational energy release was simulated by 



where L,. is the radiative flux and q = M,./M. The burning shell 
exhibited periodic flashes, all treated under hydrostatic equilibrium. 
As Fujimoto assumed steady state, the rate at which the mass of the 
WD grew was equal to the assumed accretion rate. 

Prialnik & Kovetz (1995) alleviated the assumption of steady 
state. The have shown that one must calculate at least several 
dozens of cycles before the effects of the arbitrary chosen initial 
conditions decays away. In some cases even 1001 cycles were fol- 
lowed (Epelstain et al. 2007). The calculations were extended by 
Yaron et al. (2005) with basically the same result. In particular, it 
was found that WDs accreting at rates ih > lO^^Mo/yr grow in 
mass. 

In addition to the above theoretical arguments, interesting ob- 
servational constraints should be considered as well. Some argue in 
favour of SD scenarios and some argue against. 

The SD scenarios imply that the former companion stars 
should remain after the explosion, and be observable in associ- 
ation with nearby remnants. And indeed, a G-star was claimed 
to be the former companion of Tycho Brahe's 1572 supernova 
(Ruiz-Lapuente et al. 2004). On the other hand, the opposite claim 
was also made. No former companion was found to be associated 
with SNR 0509-67.5, down to M,, = -1-8.4 (Schaefer & Pagnotta 
2012). In this respect, one should also mention the lack of any ra- 
dio detection in SN 201 Ife, which took place only 8 Mpc away. 
This was used, under several theoretical assumptions, to place an 
upper limit of order 10"*Mo/yr, on the companion wind prior to 
the supernova explosion (Horesh et al. 2012). Using Chandra and 
HST archival data, SN 201 Ife could also be used to place a limit 
on the accretion spectrum. If the WD was emitting at the Edding- 
ton luminosity, then its temperature should have been within 60 
eV > kT > 10 eV (Liu et al. 2012). In contrast, Voss & Nelemans 
(2008) reported the disappearance of an X-ray source from the lo- 
cation of a Type la supernova in NGC 1404. They estimated the 
source to have been emitting 2 - 3 x 10^^' erg/sec in X-rays. 

In terms of statistics, high rate accretion implies that SD pro- 
genitor systems should be sufficiently bright as to be easily de- 
tected. It is presently not clear at all, whether SD progenitors can be 



identified with any of the known cataclysmic binaries, the known 
symbiotic systems, or whether perhaps with the very bright su- 
persoft X-ray sources (Iben & Tutukov 1994; Gilfanov & Bogdan 
2010). 

Given the present state of uncertainty, further study of the dif- 
ferent progenitor scenarios is necessary. Here we shall address one 
particularly pressing question, which is whether systems accreting 
hydrogen at high rates can actually accumulate mass. 

We begin in §2 by describing our numerical code, and the deli- 
cate points requiring attention when carrying out a multi-cycle evo- 
lution. In §3 we describe our results. Here we will see that there are 
two types of thermonuclear runaways (TNRs), those which burn H 
and expel no matter, and larger ones burning He, which do do ex- 
pel. We will therefore discuss in §3.4 the energetics of how this is 
possible, even though the specific energy release in helium burning 
is small. We end with conclusions in §4. 



2 DETAILS OF CALCULATIONS 

Our main purpose is to investigate the long term secular be- 
haviour of high accretion rate systems. We know from the work of 
Kovetz & Prialnik (1994) that a relatively large number of cycles 
is required for the detection of the small secular effects. For this 
reason, we ensured that the code can simulate thousands of cycles. 
Even so, the calculation was carried out without compromising the 
physics. 

The initial models included WDs with masses of 1, 1.25, 1.35 
and 1.4Mo, despite the findings of Weidemann (2000) that this is 
too high and may be very rare. The models are of bare C/O WDs, 
namely, there is no atmosphere composed of lighter elements, such 
that the C/O composition extends to the surface. The initial core 
temperature is between lO'K and 5 x lO'K. However, as we will 
show later, the final results are insensitive to the core temperature. 
In addition, the models assume no rotation of the WD, and that the 
dissipation of the rotational energy at the boundary layer is com- 
plete such that the entire rotational energy is radiated away. These 
conditions are the most favourable for obtaining highly degenerate 
matter at the surface. 

All models considered accrete solar composition material at a 
rate of IO^^Mq/^t. According to Nomoto (1982), this accretion rate 
is large enough to puff the accreting system into a red-giant-like- 
star burning at steady state, with a slow evolution on a timescale 
determined by the nuclear burning shell. 

The cyclic phenomenon demands very fine division during 
mass loss and peak burning, as a coarse discretisation was found to 
introduce unacceptable errors. Therefore, the minimal mass shell 
was chosen to be between 10"' - lO^'M©, depending on the mass 
of the WD. At the peak, over 8000 mass shells were integrated. 

The code we use is based on Prialnik et al. (1981); Prialnik 
(1986). It includes diffusion, convection, and all relevant nuclear 
reactions up to A=40, which consider the intermediate electrostatic 
screening correction. 

To a large extent, mass loss following a TNR can occur 
through two types of mechanisms. The first type is through an "ex- 
plosion". Namely, a large release of energy can be used to eject 
mass until the star settles into a dynamic steady state. However, 
after a few dynamical times have elapsed, a steady state wind can 
be driven through a continuum driven wind (Bath & Shaviv 1976). 
The numerical code should be able to describe both types of mass 
loss. This can be done by employing a general mass loss condition. 

The ejection of a mass element occurs, irrespective of the 
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mechanism, if its velocity is larger than both the speed of sound 
and the escape velocity at the radius from which the mass is re- 
moved (Prialnik et al. 1978, 1979). This general rule leads to very 
extended radii and short time steps before the conditions are satis- 
fied. Kato (1983) was the first to fit a steady state wind in a model 
with two parameters which were the mass loss rate and the mass 
of the envelope. The envelope was assumed to be in place as initial 
conditions. Kato found that steady mass-loss solutions can be found 
only if the luminosity at the bottom of the envelope is greater then 
the minimum value of the local Eddington luminosity throughout 
the envelope. 

Here we follow Kovetz (1998), and check at every time step 
whether an optically thick wind can exist, and remove mass ac- 
cordingly. We stress that in no case, did the removal of the mass 
was found to induce additional mass removal, namely, no rarefac- 
tion wave was created leading to endless mass loss. The possibility 
of an optically thin wind was not considered in this calculation. 

We quench the accretion during the eruption, because we ex- 
pect the accretion disk to be destroyed under the dynamic and high 
luminosity conditions. Thus, mass accretion is switched off when 
the nuclear luminosity is larger than 10 times the total luminosity, 
and switched on again when the luminosity decreases below the 
threshold' . For this reason, the long term average accretion is typi- 
cally a third of the peak instantaneous accretion rate of lO^^Mo/vr. 



3 RESULTS 

Figure 1 depicts the initial behaviour of the I.OMq WD model, 
which settles into steady cycles of small outbursts without eject- 
ing any mass. Unlike accretion at low rates, giving rise to classical 
or recurrent novae (Rakavy & Shaviv 1968), the high rate accretion 
keeps the accreted layers barely degenerate. As a consequence, the 
TNR is modest and involves a small shell. 

Figure 2 demonstrates that the cyclic behaviour continues for 
many cycles, with a period of about 10 yr. There is very little secu- 
lar change from one cycle to the next. In each cycle, hydrogen un- 
dergoes a small runaway in which the luminosity increases by over 
a magnitude. The flashes are characterised by a very sharp rise and 
a slower decline. A refined picture is shown in panel B of fig. 2, 
demonstrating how the system approaches the Eddington luminos- 
ity for a couple of years. The system stays in the minimum state 
only a very short fraction of the total period, until thermal ignition 
starts a new cycle. 

Because the emitted luminosity during the flash is slightly less 
than the Eddington luminosity, no mass loss develops. This is be- 
cause the luminosity was used to expand the thin hydrogen shell 
(see fig. 3). 

The burnt hydrogen accumulates as helium. The abundances 
profile after 3000 cycles of small outbursts are shown in fig. 4. The 
underlying CO WD is present at small radii, below ~ 7 x lO^^Rg. 
It is surrounded by a thick layer of helium with smaller amounts of 
C and O. The freshly accreted material having solar composition 
is located at the outer radii. It is important to note that the burning 
is almost complete. At this phase the system appears as a perfect 
recurrent object which accretes the ashes onto the original WD. 



' Using other conditions, such as switching off the accretion only when 
mass loss is present gave rise to minor changes in the time scale. 
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Figure 1. The first 10 outbursts for a I.OMq WD with peak accretion rate 
of lO^^Movr"' - only the first outburst ejected mass. 
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Figure 2. The secular sequence of flashes for the model depicted in fig. I , 
3 X lO** years after the onset of accretion. Panel B shows one cycle in detail. 
Note the high similarity from cycle to cycle. 



3.1 Secular changes in the accreting star 

The initial profile which we employ is obtained by evolving a poly- 
tropic WD, until a desired central temperature is reached. Once 
accretion cycles begin, the small hydrogen flashes give rise to a 
temperature wave which propagates into the helium layer and the 
outer edge of the CO WD. This wave increases the temperature at 
the base of the accreted helium envelope to nearly 1.2 x lO^A^ after 
4000 cycles. Fig. 5 provides snapshots of the temperature profile. 
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Figure 3. The maximum radius reached during outbursts, for the model 

depicted in fig. 1. The small variations in the maximum radius are due to Figure 6. The temperature profile at the core, for the model depicted in 

the finite time step taken. gg j Note that the heating is due to the adiabatic contraction of the WD 

associated with the additional accreted mass. 
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Figure 4. The abundances after 3000 cycles of small outbursts, for the 
model depicted in fig. 1 . 

The large mass of the accreted envelope, nearly 0.0 IMg, af- 
fects the structure of the WD as well. As a consequence, the 
WD contracts and the inner temperature rises from 6.0 x 10' A" to 
6.17 X 10' A"— a 3% change (see fig. 6). 

It is interesting to compare the temperature profiles of the 
WD for difi'erent accretion rates. We therefore, calculated using 
the same initial WD structure a model having an accretion rate 
of lO^'MoV"'. We follow the calculation for 3000 cycles. Note 
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Figure 5. Temperature profiles at different cycles, for the model depicted 
in fig. 1. The slow diffusion of heat into the WD is due to the temperature 
inversion which inhibits convection, and the fast time scales associated with 
the high accretion rate which limits the energy lost outwards. 
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Figure 7. Comparison between the temperature profiles after 3000 outbursts 
of a 1 Mq WD, accreting at a rate of lO^'Mo/.yr (dashed) and at a rate of 
IQ-^M^Iyr (solid). 



that in both case, this corresponds to a total accreted mass of order 
IQ-^Mq. See fig. 7. 

For the low accretion rate, where the accreted material is suf- 
ficiently degenerate, the outbursts are strong enough to eject all the 
accreted material, as well as some of the eroded WD. The interval 
between outbursts (3 x \(fyr) is long, such that the WD can cool 
down. The temperature profile of the WD accreting is then hardly 
affected by the outburst itself. The only evidence for the existence 
of the previous outbursts is a small hump in temperature profile. 

The outcome of the high accretion rate is completely dif- 
ferent. Since the material is insufliiciently degenerate, the pres- 
sure at the base of the envelope prior to the outburst, P„i, ~ 
GM,„iim„,y/4KR^^,j s; 8 - 9 X 10"dyn cm^^, is too low to allow mass 
ejection. This result is in agreement with the analytical calcula- 
tion of Fujimoto (1982) and MacDonald (1983) that predicted mild 
outbursts occur for P„ii ;£ lO'^dyn cm^^. Because the outburst is 
weak, matter cannot escape the gravitational field of the WD. Con- 
sequently, a heat wave can form and penetrate into the base of the 
accreted envelope. This heat wave will serve as the trigger for the 
giant outburst described in the next section. 

To summarise, the positive temperature gradient in the high 
accretion rate inhibits heat flow into the core. The degree of heating 
per cycle is negligible, but accumulates during thousands of cycles. 
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Table 1. The total accumulated mass as function of WD mass for M 
10"*Mo>'r"' before completely ejected by the final outburst 
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Figure 8. The nuclear luminosity during the final outbursts. 



3.2 The final outburst 

The surprise came after 4153 cycles when the accreted mass 
reached IO^^Mq. The maximum temperature in the helium layer 
reached 1.2 x lO^K, allowing helium ignition in a strong flash (see 
figure 8). 

The peak nuclear energy generation during the He TNR 
reached lO'^erg/sec and the maximum temperature reached a peak 
of 4x lO'K. 

The very strong He flash caused a temperature wave which 
propagated to the hydrogen burning shell located above the pure 
helium layer. The wave raised the temperature, and accelerated the 
burning accordingly. This resulted in the ejection of about 1/3 of 
the helium layer accreted. This phenomenon repeated two more 
times with the net outcome that the entire accreted helium layer 
was ejected. 

These results should be compared with those of Yaron et al. 
(2005), who calculated a similar model. For m„cc = lO^^Mo/jr, 
they have found that the WD erodes. For = lO^^Mo/vr, 

the white dwarf gained mass through outbursts, while for ot„cc = 
10"*Mo/yr it gained mass without outbursts. However, they only 
followed a limited number of cycles, which did not allow for He 
ignition. 

We repeated the calculation for colder WDs, having lower core 
temperatures. We found minor changes. For example, the time be- 
tween outbursts increased from 8 yr for a 6 x lO'K WD to 9.5 years 
when the central temperature was lO'K. He ignition took place af- 
ter 4900 accretion cycles, corresponding to an accreted mass of 
0.0103Mo. 

We also repeated the calculation for WDs having different 
masses. The time between outbursts varied from 8 years for a IM© 
to three months for the 1.35Mo. But the main results remain the 
same. The results are summarised in Table 1. In all three calcula- 
tions, the last stage was nearly the same — once the pressure at the 
base of the He envelope became sufficiently high, the He ignited 
and the entire accreted envelope was ejected. 

Unlike lower mass WDs, the 1.4Mq WD behaved differently. 
During the initial hundred or so cycles, it followed a similar be- 
haviour as lower mass WDs — small outbursts took place with a rep- 
etition period of ~ 18 days. However, after roughly 100 cycles, the 
WD expand and became a red giant. The initial cyclic behaviour we 
find is consistent with Kovetz & Prialnik (1994), who investigated 
the effect of accretion onto a 1 .4Mq CO WD at different accretion 
rates, between 10"'" to lO'''Mo/yr. 



WD Mass [Mq] 


Final accreted mass [Mq] 


1.0 


9.9 X 10"^^ 


1.25 


7.5 X 10-'' 


1.35 


1.1 X 10-'' 



Table 2. Ejecta Abundances. 



Element 


Solar 


Maximal mass 


H 


0.7 


6.3(-5 ) 


"He 


0.28 


0.71 


12q 


3.9(-3) 


1.25(-1) 


14n 


1.0(-3) 


5.4(-3) 


160 


9.4(-3) 


5.7(-2) 


2«Ne 


0.0 


4.3(-2) 


"Na 


0.0 


1.5(-4) 


^^Mg 


0.0 


5.9(-2) 


27 Al 


0.0 


1.2(-5) 


28Si 


0.0 


3.8(-4) 



3.3 Abundances 

The abundances in the mass ejected in the final three outbursts are 
given in Table 2. We find the following. 

(i) There is very little hydrogen ejected. This is because most of 
the hydrogen was converted into helium during the small outbursts. 
The small amount which is ejected is from the outer most accreted 
layer which did not have time to burn. 

(ii) Helium abundance increased to about 0.7. This is more than 
twice the solar composition. However, it is less than the initial X+ Y, 
because of the triple a reaction. 

(iii) The Z increased by more than an order of magnitude, from 
0.02 to nearly 0.29. The interesting point, however, is that the iso- 
topic ratios are different from those in the standard solar composi- 
tion (or different for standard nova ejecta). 



3.4 Understanding the energetics 

A simple analytic approximation for the radius of a zero tempera- 
ture WD was given by Verbunt & Rappaport (1988), it is 



0.0114- 



((M„,rf/Mrf,)-2/3 _ (M„rf/Mrf,)2/3)"' 
(l + 3.5(M„,rf/M„)"" + {M„JM„y'J" 



(1) 



where M,, = 0.00057Mo is a numerical constant. Using this ex- 
pression, we can calculate the critical mass M„„ above which the 
binding energy of an element is larger than the specific energy e 
released in nuclear burning. Namely, we require that 

GM„i, 



(2) 



where ii/j,,,^ is the nuclear binding energy per unit mass. The results 
are given in table 3. 

According to the standard lore, hydrogen burning can there- 
fore eject mass from WDs having masses of up to 1.436Mo. On the 
other hand, helium burning can eject mass provided the mass of the 



© 2010 RAS, MNRAS 000, ??-?? 



6 Idan et al. 



Table 3. The critical mass Mcru above which the binding energy is greater 
than the energy released in nuclear burning. 



Reaction 


Energy release [MeV] 


Q [erg/gr] 


Merit 


Ap -> He 


26.47 


6.3 X 10'** 


1.436 


3He ^ C 


7.276 


5.8 X 10''' 


1.30 


C + C ^ Mg 


7.075 


2.8 X 10'^ 


1.29 



WD is less then I.BMq. Nevertheless, this is the opposite of what 
the present models exhibit. 

In all our 10"'Mq/>t models, with WD masses up to 1.35Mq, 
the hydrogen flashes did not expel any material even though their 
mass was below the critical one. On the other hand, all these models 
did expel the helium accumulated from the hydrogen burning, even 
though the most massive of which were already above the critical 
WD mass for helium ejection. 

To understand how hydrogen burning could result with no ex- 
pulsion, but He can, let us consider the energetics. 

When hydrogen is initially accreted on the bare WD, at an 
average rate m, it releases its binding energy (GM/Ry,/o)ih. In ad- 
dition, as it burns, it releases riien, where £u is the specific nuclear 
energy released in hydrogen burning. On the other hand, the max- 
imal luminosity with which the system can lose energy is tjLejj, 
with being some factor smaller than unity. This factor depends on 
the actual luminosity of the expanded state, and its duty cycle. For 
the IMo WD accreting at lO^^Me/jf-, it is ^ ~ 0.5. 

Comparing these two energies, gives m(GM/RwD + €u) > 
TjLEdd- In other words, every H burning cycle releases more energy 
than can be radiated away. As a consequence, the envelope must ex- 
pand to some effective radius, for which the specific binding energy 
satisfies 



IGM 



(3) 



Moreover, the gas in the expanded state is hot, and therefore its 
specific internal energy Uj„, is important, and can be as much as 
half the gravitational binding energy. 

Once helium is ignited, the specific binding energy of the 
puffed up envelope is therefore smaller than the specific binding 
energy it would have had if the helium in the envelope was degen- 
erate and located at r a Rwd- Thus, when helium ignition takes 
place, the specific nuclear energy released is sufficient to expel the 
envelope. Namely, 



GM 

Rei 



+ Him + ^He > 0. 



(4) 



Clearly then, the main reason why the helium burning can expel 
matter is because the high accretion rate forced the accreted mat- 
ter to expand and heat such that it had a smaller specific binding 
energy. One can recast this difl'erently. If we look at the average 
binding energy, -GM/R^ff + ;(,„, from eq. 3, and plug it in eq. 4, we 
find that the condition for expulsion is: 



(5) 



Thus, the high accretion rate, the hydrogen burning, and the limited 
luminosity combine to keep the envelope bloated enough to allow 
expulsion. However, an accretion rate which is too high would force 
the system to become red-giant-like, and the envelope would be 
expelled. 
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Figure 9. The behaviour of the convection zone duiing a hydrogen flash. 
The broken line represents the temporal evolution of the C-O/He boundary. 
The dotted line describes the behaviour of the He/H boundary. The upper 
line is the radius of the outermost shell. 
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Figure 10. The profile of the specific entropy a at the begin- 
ning of the 3000''' cycle. The entropy is given in units of sq = 
7T-(MalNA)- (GincfK-^. 



Note also that eq. 3 should include two additional components. 
First, some of the heat in the envelope should diffuse into the WD. 
This implies that there is an additional thermal energy term which 
should in principle appear. However, when this energy was artifi- 
cially omitted, the final results hardly changed. 

Second, an additional energy gain by the system is associated 
with the increased binding energy of the underlying WD — the lat- 
ter contacts and adiabatically heats. However, since the heat cannot 
diffuse out over the n accretion cycles, it is stored and released upon 
the final ejection. In other words, this energy component is unim- 
portant. 



3.5 Few more results 

The history of the convective zone is depicted in fig. 9 along with 
the borders between the hydrogen and helium layers, and between 
the core and the helium layer. In fig. 10, the profile of the specific 
entropy is shown at the beginning of the 3000'* cycle (see fig. 9). 
Note the the entropy profile is given only for the inner region where 
the convection did not penetrate. As is obvious, there are two en- 
tropy barriers which prevent the convection from penetrating any 
deeper and mix the C and O extensively over the envelope. The 
resulting abundances are shown in fig. 4. 
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Table 4. Characteristics of the Nova Envelope for Accretion Rate of 10 M^yr 





Core Temp. 


Accreted Mass 


Ejected Mass 


Residual Mass 


Recurrence 


[Mo] 


[\(fK] 


Per Cycle [Mq] 


Per Cycle [Mq] 


Per Cycle [Mq] 


Time [yr] 


1.00 


60 


9 X 10"* 


8 X 10"* 


1 X 10"* 


100 


1.25 


50 


1.5 X 10"* 


1.2 X 10"* 


3 X 10"' 


16 


1.35 


50 


3.5 X 10"' 


2.6 X 10"' 


9 X 10"* 


3.7 


1.40 


60 


5 X 10"* 


4 X 10"** 


1 X 10"* 


0.5 



3.6 Accretion rate of 10 'Mq>t ' 

According to Yaron et al. (2005), WDs should erode for accretion 
rates lower than about 10"*Mo,yr"'. However, as we have seen 
above, the option of achieving the Chandrasekhar limit is also ex- 
cluded for an accretion rate of 10"*Moyr"', because of the strong 
helium flashes that take place after several thousand hydrogen igni- 
tion cycles. We have therefore decided to include the simulation of 
WDs accreting at an intermediate rate of 10"'Mo,yr"'. We found no 
helium flashes. Moreover, since the conditions do not satisfy eq. 5, 
even if there were, there would have been no expulsion. The results 
of these models are summarised in table 4. 

Our results confirm those of Yaron et al. (2005). Namely, the 
WD experiences a net increase in mass when accreting at a rate of 
10"'Moyr"'. However, the residual mass fraction left on the WD 
is relatively small, of order 10-25%. Although attempts to identify 
such a system with known observed ones failed, such systems if 
existed could in principle approach the Chandrasekhar mass. 



4 CONCLUSIONS 

High rate accretion of solar composition material onto WDs pre- 
vents degeneracy. According to general wisdom, this should give 
rise to quiet steady state burning in which the helium ash from the 
hydrogen burning in the envelope should accrete onto the core, and 
secularly increase its mass. 

Since eruptions are expected to disrupt the accretion process, 
our first conclusion is that instead of a steady state process, the 
accretion and burning take place in eruptive cycles. However, this 
cyclic behaviour is different from classical nova eruptions. We find 
that the accreted material is only moderately degenerate before 
it ignites. And once it does, it burns completely. Therefore, any 
ejecta will contain no detectable hydrogen. Clearly then, the non- 
existence of hydrogen in observed ejecta does not immediately im- 
ply that the accreted matter did not contain any hydrogen as well. 
Since the lack of observable hydrogen in ejecta is generally con- 
sidered as one of the strongest evidence suggesting the accretion of 
pure helium, such claims should be considered more cautiously. 

The last, but most important conclusion we have reached is 
that high accretion rates, of 10"*Mo, do not lead to the a secu- 
lar increase in the mass of the white dwarf, because giant helium 
eruptions take place after a sufficient amount of helium accumu- 
lates. These eruptions expel all the accreted layer and even a small 
amount of the underlying WD. 

We find that helium flash explosions take place for WDs with 
masses up to 1.35Mo. From eq. 5, however, there should be no up- 
per limit above which helium cannot expel the material — although 
more compact systems would require more energy to expel material 
out of them, they also release more energy during the accretion. 

The above accretion rate is very high and can be maintained by 



only few systems. The results for a lower accretion rate, of 1O"'M0 
is qualitatively different. There are no giant helium flashes, and 
there is a net increase of the WD mass. However, because the mass 
accretion efficiency is small, it requires the donor to fine tune the 
mass transfer over a large integrated total mass. It therefore seems 
unlikely that large increases in the WD mass are possible. 

The above conclusions depend on the energy release by the 
hydrogen burning. Thus, systems in which helium is accreted will 
not have their envelope bloated, and any helium eruption will not 
lead to large mass loss. This possibility is presently under investi- 
gation. 

Another important point that is the subject of further inves- 
tigation is the effects that will arise once super-Eddington condi- 
tions are allowed to develop (Shaviv 2001). This should modify the 
above conclusion for two main reasons. On one hand, the super- 
Eddington luminosities drive mass loss through continuum driven 
winds. On the other hand, a higher luminosity will inhibit the at- 
mosphere from bloating, such that the helium flashes will not nec- 
essarily be able to drive the envelope away. Thus, it is not a priori 
clear how the above picture will be modified. 
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